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Abstract. Many studies on animal and human movement patterns report the 
existence of scaling laws and power-law distributions. Whereas a number of random 
walk models have been proposed to explain observations, in many situations individuals 
actually rely on mental maps to explore strongly heterogeneous environments. In this 
work we study a model of a deterministic walker, visiting sites randomly distributed 
on the plane and with varying weight or attractiveness. At each step, the walker 
minimizes a function that depends on the distance to the next unvisited target (cost) 
and on the weight of that target (gain). If the target weight distribution is a power-law, 
p{k) ~ , in some range of the exponent /3, the foraging medium induces movements 
that are similar to Levy flights and are characterized by non-trivial exponents. We 
explore variations of the choice rule in order to test the robustness of the model and 
argue that the addition of noise has a limited impact on the dynamics in strongly 
disordered media. 
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1. Introduction 

Deterministic walkers, as opposed to random ones, follow non-stochastic rules of motion. 
They can be used as an approximation to describe mobile agents that (i) have some 
information on the medium they explore (composed, for instance, of food patches, 
facilities or cities, depending on the context); and (ii) use this information to optimize 
their gain (such as food intake or the possibility of realizing a particular activity). If the 
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medium is disordered and heterogeneous, deterministic rules can produce fairly erratic 
and complex trajectories, susceptible of being studied statistically as random walks. 

Unlike Brownian particles, humans and other animals keep memory of their past 
activity. For example, a class of deterministic walks assumes that each decision depends 
on the previous steps the walker has performed. This kind of processes can be applied to 
a wide range of hard-to-solve problems such as the travelling salesman or the travelling 
tourist [Ill2]- Yet, despite their potential importance in a variety of problems in physics, 
mathematics, computer science, economics and biology, the properties of deterministic 
walks remain largely unexplored. Deterministic diffusion has been investigated in simple 
random media [31 H]. Similar processes can be used to implement classification schemes 
[5l [6] and for pattern recognition [7]. In biology, though not widely used, deterministic 
walks offer an appealing framework to model the movement of foraging animals. Indeed, 
many organisms rely on mental maps to navigate their environment in a non-random 
way [8]. Primates and other animals dispose of sophisticated cognitive skills [9], such 
as the ability to remember autobiographical events that happened at some particular 
places and times (episodic memory). 

The presence of scaling laws and power-law distributions in animal movement 
patterns is an ecological problem that has attracted considerable attention in the last 
decade [10] . Empirical evidence of Levy- walks has been reported, among others, for 
micro-organisms [H], bumblebees [12], spider monkeys [I3], marine predators |14j . 
Similar patterns have also been observed in humans, as shown by studies on hunter- 
gatherers [15], bank notes [16] or fishermen boats Human individuals also display 
multiple scale displacements well described by truncated power-laws [18j. Most of the 
theoretical interpretations proposed to explain observations rely on the random walk 
hypothesis or some variant thereof: for example, the diffusion of bank notes is actually 
well described by a continuous time random walk model [T6] . 

A common hypothesis is that the ubiquitous presence of widely fluctuating 
characteristic scales in animal movement could optimize the success of search processes 
of randomly distributed resources [12]. According to this perspective, animals would 
execute a particular optimal foraging strategy, usually a Markovian stochastic process 
(see [121 [ini EO] for examples), to find preys that are uniformly distributed and only 
detectable at short distance. These strategies are susceptible of being learned along an 
individual life experience [T^, or transmitted through generations as an evolutionary 
trait [21]. 

An alternative hypothesis, namely, that scale-free (or more generally, multiple scale) 
movement patterns may be an outcome of the direct interaction of the forager with a 
complex environment, has been much less explored quantitatively. In this approach, very 
simple movement rules, sometimes involving memory, may produce complex searching 
patterns. For instance, a simple response to passive scent concentrations in turbulent 
structures can account for scale-free trajectories of flying insects [22]. Levy patterns 
similar to those observed in jackals can also arise from avoidance of conspecific odour 
trails, as modelled by systems of self-avoiding walks [23]. Levy- like distributions of 
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the distances between detected preys can also emerge in the trajectories of a predator 
following, deterministically, chemotaxis gradients produced by immobile and randomly 
located preys [21]. 

In this article, we study a model describing the movement of a walker that forages in 
a strongly heterogeneous environment using mental maps. It is assumed that the forager 
has a perfect knowledge of the available resources (for instance, fruiting trees) and that 
it chooses deterministically at each step the "best" (in terms of an efficiency function) 
unvisited food patch in the whole system. Therefore, decisions are not taken on the 
basis of the knowledge of the immediate surroundings only. We originally proposed this 
model in the context of a study on the ranging patterns of spider monkeys in tropical 
forests [25l[26], but we expect it to have a much wider application range. In that work, 
we showed that scale-free distributions of step lengths can emerge if the resources are 
sufficiently diverse in size. Here we give a detailed analysis of the model, we show that 
the behavior of the walker is robust to the addition of noise and we study how the results 
are affected by changes in the deterministic decision rules. 

2. The model 

Consider a two-dimensional square domain of unit area, randomly and uniformly filled 
with point-like targets of fixed positions. The targets model, for instance, fruiting 
trees in an ecological context, or cities in a social one. The targets may also represent 
network communities, nodes for information retrieval or markets for salesmen. Each 
target i is characterized by, say, a food content or a size, ki, which are independent 
random integers drawn from a distribution p{k). (Integers are considered here for 
computational convenience, without restricting generality.) The distributions of tree 
sizes in tropical and template forests are very broad and in many cases well described 
by inverse power-laws with exponents ranging from 1.5 to 4 [271 [28]. Across species, the 
size of a tree is roughly proportional to its fruit mass [29]. We will therefore consider 
target size distributions of power-law form: 

p(k) = Ck-^, with = 1, 2, kraa.] > 1, (l) 

C being the normalization constant; and kmax is a cutoff size, that can be taken infinite 
(see section [6]) but will be fixed to 10'^ in the simulations below, except when otherwise 
indicated. In the social context, ^ with /3 = 2 is equivalent to assuming that city sizes 
are distributed according to a Zipf's law [30] . 

In this domain, a forager (or traveler) is initially located on a target. We assume 
that the forager has a perfect knowledge of the size and position of every target in 
the system. In a human context, when the travelers are, say, salesmen or tourists this 
assumption is natural. On the other hand, in a biological context, this assumption is 
based on observations that many animals (primates [HI [9], birds [31]) use local and global 
cues during searching. Field studies show that Clark's nutcracker (among other species 
of small birds) can store 22,000 to 33,000 nuts in up to 2700 locations over an area larger 
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than 100 square miles and remember where almost 70% of them where placed [32] • One 
could also relax the perfect knowledge hypothesis by assuming that the forager knows 
a random subset of Ns targets, or the Ns largest targets {1 Ns < N). These known 
targets would still be described by the distribution ([1]) with exponent — /3. 

At each time step (t — *• t + 1) the following rules are iterated: 

(i) the walker located at target i will next visit the target j such that the step 
efficiency Eij is maximum among all allowed j 7^ i in the system, 

(a) previously visited targets are not revisited (which corresponds to destructive 
foraging in the ecological context). 

Presently, we will consider two forms for the efficiency function Eij: 

Eij = kj/lij ("quotient" rule) (2) 
Eij = kj — hj/lo ("difference" rule), (3) 

where kj is the size of target j, kj the distance separating targets i and j, and Iq = 1/ \fN 
the typical distance between nearest neighbour targets. Hence, in both rules (I2])-(I3]), 
the forager tends to visit the largest possible target accessible in a minimum travelled 
distance from its current position, ignoring targets that are either too small or too far 
away. Rule (I3l), which has not been studied before in [26], represents the net energy 
balance of a sojourn, assuming that the energy expenditure is proportional to the 
distance travelled. Note that in this case we do not impose the condition Eij > 0: 
some steps may be such that energy expenditure exceeds energy gain (see section [5] 
below) . 

3. Step lengths statistics 

In [26], a mean field, single-step solution of the model for the quotient rule ([2]) was 
obtained in any spatial dimension d. There, it was shown that the heterogeneity of the 
medium can induce power-law distributions P{1) of sojourn lengths (distances between 
the origin and the best target in that case): In (i = 2, if /? > /?c = 3 in the target size 
distribution ([1]), then 

P{1) oc with aMF = (3-2. (4) 

Therefore displacements can be considered as Levy-like (1 < a < 3) if 3 < /3 < 5. If 
/3 < 3 and kmax = 00, the density of very big targets in the medium is so high that P{1) 
is no longer scale-free and involves a characteristic size diverging with [26] . 

Although the mean-field description captures the gross behaviour of the model, 
it fails to explain the numerical results quantitatively (see below). By imposing the 
condition that the walker cannot return to a previously visited site, memory effects 
produce non-trivial long range correlations between steps as the process develops. In 
[25] , 2d simulations with the quotient rule ([2]) indeed exhibit scaling-laws, but different 
exponents are found: 

P{1) oc with a ~ /? - 1(^ omf), for 3 < /? < 4. (5) 
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Figure 1. Left panels: Examples of deterministic trajectories for 3 values of 
the target size distribution exponent: f3 ~ 1.5 (superabundant medium), 3 (diverse 
medium) and 6 (homogeneous medium). When f3 — 3, the sojourn length fluctuates 
widely in both cases, similarly to a random Levy flight. {N = 10^, number of steps 
per walk= 10^.) Right panels: Step length distributions. Near-perfect power-law 
distributions, P{1) ~ are obtained in the range 3 < /3 < 4. The value of the 
scaling exponent a,t P — 3.1 is a ~ 2.1 while the value at /3 — 4.0 is a ~ 3.0 

We discuss in more detail the numerical results below. 

Several trajectories are represented in figure [T]-left, with the corresponding step 
length distribution in figure [T]-right. With the "difference" rule ([3]), three qualitative 
regimes can be identified as j3 is varied. In the first one (superabundant regime), 
corresponding to /? < 3, the trajectories are dominated by large sojourn lengths. If 
/3 is small enough, there are many attractive (large) targets, so that the forager often 
travels large distances to reach them. Very large targets are relatively distant from 
each other, but the large values of their size k make these long sojourns "affordable". 
The situation changes dramatically in the interval 3.1 < j3 < 4.0 (diverse regime), 
where very big targets exist but are too scarce to produce frequent long travels. As a 
result, there are a few large sojourn lengths alternated with many shorter ones at all 
scales. In this regime, the distribution of step lengths is a power-law with an exponent 
a approximately given by relation ([5]). Trajectories in this regime closely resemble Levy 
flights. Finally, when (5 > 4 (homogeneous regime), target size fluctuations are too 
small and the forager performs mostly short steps to nearby targets, similar to a (self 
avoiding) random walker. 

The three regimes discussed above, in particular the scaling regime with a ~ /3 — 1, 
see ([5]), are also observed when the walker follows the "quotient" rule ([2]). In the 
superabundant regime, nevertheless, the steps are on average shorter than for the 
difference rule (see flgure [TJ-right). 
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Figure 2. Length fluctuation ratio, as a function of the resource exponent 

/3. In both the k — I/Iq and k/l rules, this quantity is maximum at f3 ~ 3.1. In all 
cases, N = 10^ and kmax = 10'^. Data was obtained from averaging over 10 different 
disordered media. In each realization, a walk starts near the center of the square 
domain and visits 10^ targets, a number still <C N. 
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Figure 3. Iq'^ /{l)^ as a function of /3 showing diverging steps below /3c. kn 
other system settings as in figure [2] 
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Further insight into the statistical properties of the walks can be gained by studying 
step length fluctuations. Figure [2] displays the fluctuation ratio, deflned as (/^)/(/)^, 
versus the resource exponent f3. Fluctuations are maximum at /?c ~ 3 (independently 
of kmax) for both decision rules, conflrming the results obtained in [25] for rule ([2]) only. 
Therefore, a medium with (3 = (3c produces, in some sense, the most scale-free walks, 
that are characterized by a ~ 2. As shown in flgure[3], f3c also corresponds approximately 
to the onset of emergence of very large steps: (/q/ (0)^ becomes very small below jSc- 
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Figure 4. Temporal evolution of the average size of the i visited target, in the k — I 
rule and with (3 = 2.25, N = 10^, kmax — lO"^- Averages are taken over 10 independent 
runs. 
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Figure 5. Average time derivative (fc^+i — ki)i performed over the time domain 
i = 1,...,10'' (for the rules k — and k/l) as a function of (3. Each point is the 
average of ten replicas. The quantity ki is defined as the average value of r consecutive 
ki values. Notice the onset of stationarity at (i — i. 



4. Sizes of visited sites and non-stationarity 

We now analyze the size of visited targets, which is related to the intake of the forager 
in the ecological context. Figure H] shows the temporal evolution of the average size 
{ki) of the i^^ visited target, for the k — 1/Iq rule and a superabundant medium with 
j3 = 2.25. For this value of (3, it is clear that this quantity is not stationary. Three 
qualitative stages can be distinguished. The initial stage is characterized by the visit 
to large targets all of which are relatively similar in size; in a following stage, the size 
of visited targets decreases slowly with time; in the final stage (where still 99% of the 
targets have not been visited) the visited targets have small sizes. 

This behaviour is typical of the superabundant regime (/? < /3c) and can be 
understood qualitatively in terms of the length scales considered by the forager. In 
the initial regime, decisions are taken essentially in terms of the target size, the distance 
travelled playing a minor role (see section [6] for a analytical argument). In this regime, 
large targets are so big that the walker chooses them in spite of having to take long steps 
to reach them, steps that can actually be of the order of the system size. Once the large 
targets have been depleted, there are no targets left that are worth travelling distances 
of the order of the system size, but there are still plenty of targets worth travelling 
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Figure 6. Left panel: Average energy balance as a function of /? for three values of 
A in rule ([6]). Insets: Eij{t) as a function of t for [3 — 2> and 4 (A = 1 in both cases). 
Right panel: Step length distribution at fixed (3 — i and for three values of A. 



distances which are much longer than the typical inter-target distance Iq. Finally, in 
the last regime most of the valuable targets have already been depleted, and the walker 
chooses among the many small targets that are at distances of the order of Iq. 

Figure [5] shows the average decay rate of {ki). It is interesting to note that the 
process becomes essentially stationary for /5 > (3c. Therefore, in the scaling regime ([5]), 
even though the forager modifies its environment, the medium explored is statistically 
independent of time (as long as the number visited target is -C N). 



5. Energy balance and power-laws 

Some insights into the properties of the model with the "difference" rule ([3]) can be 
obtained by modifying the energy cost associated with travels. We introduce a new 
parameter, A, which is the energy spent by the forager per unit length Iq: 

Eij = kj - Alij/lo. (6) 

In general, we expect that different values of A will induce (in a same medium) changes 
in the target choices, that will translate into changes in the length of the steps and 
therefore in their distribution. 

The case A = 1 is the one analyzed so far. If A ^ 1, we expect the average energy 
balance (Eij) to be negative (despite of the optimal decisions) and the forager starves. 
When A <^ 1, on the other hand, we expect the energy gained to be much larger than 
the energy consumed in travels. The distances between successive targets should then 
not play a major role in the decisions. Figure El left panel, shows {Eij) as a function of 
P, in the stationary regime {i.e. (3 > 3), for A = 1, 5 and 0.2. We observe that the case 
A = 1 is of particular interest and also biologically relevant, in the sense that {Eij) > 
and {Eij) = 0(1) for all values of (3 > 3. These properties do not hold in the other 
two cases. For A = 5, {Eij) < for all values of (3 > 3: a forager with such energy 
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expenditure would not survive in these media. For A = 0.2, we find (Eij) is quite larger 
than 1 at /5 = 3. 

In the insets of the left panel of figure Ej we display typical time series Eij{t) for 
A = 1 and two different values of /?. We observe in both stationary process. 

Remarkably, temporal fluctuations in the case (3 = 3 contain a few bursts of size much 
larger than the average, a feature reminiscent of a critical phenomena. 

The right panel of figure [6] shows the step length distribution P{1) at /3 = 3, for the 
three values of A considered above. Only the case A = 1 clearly exhibits a power-law 
(oc /~^). We conclude that power-laws emerge in this model from the interplay between 
the energy cost due to travels (that involves a space variable, Uj) and the energy gain 
{kj, which is not a space variable), when both terms are on average of the same order. 

6. The effect of noise in the choice of a target 

We analyze in this section a probabilistic extension of the deterministic decision rule, 
based on the "quotient" efficiency ([2]). The approach presents some similarities with 
that of [331 El] where a stochastic version of a deterministic walk visiting point-like 
targets was studied. 

In principle, different approaches can be followed to introduce noise in the model. 
In one such approach, studied in [25], the forager still seeks to optimism its movements 
but has an imperfect knowledge of the medium {e.g, the values of tree positions and 
food contents in the mental map differ from the real ones by some random fluctuating 
amounts); hence, the efficiency of a step is not evaluated exactly and the trajectory 
can depart from the ideal deterministic one. For this case, it was shown that in the 
Levy regime these forager mistakes have little effect on the step-length distributions 
[25] . Here we focus on a different approach for introducing noise into the system. In the 
present approach the forager still has a perfect knowledge of the medium and correctly 
evaluates step efficiencies, but it may decide not to visit the best target. This mimics 
the free-will, so to say, of the agent, allowing for non-optimal choices. 

For simplicity, let us consider a one-dimensional lattice with N regularly spaced 
targets at coordinates i = 1,2,. ..,N and a walker initially located at coordinate 0. 
Similarly to ([1]), target weights are quenched and distributed according to 

Pik) = i(3-l)k-^, (3>1, (7) 

with 1 < A; < oo now a continuous variable. We now assume that the walker chooses 
a target i with a probabihty pi that depends on the efficiency Ei = ki/i of the step. 
We assume pi oc f{Ei), with f{x) an increasing function of x. Hence, high efficiency 
targets have higher probabilities of being chosen, but the process is not deterministic. 
For further simplicity, we study the single step version of this problem. A somewhat 
similar diffusion problem has been studied in the weak disorder limit in [35]. Similar 
ideas were also applied to the study of a glass transition in a kinetic trap model with 
hopping rates depending on inter-trap distances [331 El] ■ 
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Let us denote rj = f{Ei), hence, Pi = ri/^.Vj. In this section we calculate the 
overlap defined as the average probability that two independent walkers located at 
the origin in the same medium choose the same target: 

N 

Yi = {y2), with 2/2 = XI P'' (8) 

1=1 

where (.) denotes the average over disorder configurations {ki}. A purely deterministic 
process corresponds to I2 = 1, whereas a purely random choice leads to 12 ~ 1/^- 

As shown in two examples below, this problem is similar to the random energy 
model of spin glasses [36]. Whereas I2 — > as ^ cxd in the disorder-free homogeneous 
case [ki = constant), there is a 0(1), or even unit probability that the two walkers choose 
the same site in strongly heterogeneous media (/? < 2, here). 

First, we note that 

sr^N 2 ^ roo 

2/2 = ^'=' 2 =J2 dttrje-'^-^e-'^^^^'^ (9) 



i=l 



where the identity A = dt texp{—tA)/T{2) has been used. Since the weights are 
independently distributed among targets, then: 

Y, = Y, dtt{rje-'^^)\[{e-'^^). (10) 
We analyze below two simple choices for the preference function f{x). 
6.1. Exponential preference 

If f{x) = exp{x/T), where T is a constant, then rj = exp{ki/iT) and the deterministic 
case is recovered as T — > (zero "temperature"). At finite T, we decompose ( ITO!) into 
Y2 = Jq dt{..) + dt{..) = + , with a a constant <^ 1. To evaluate the integrand 
of the first integral, we first approximate exp(— trj) to 1 if < 1/t, and to otherwise. 
Therefore, 

/-iT Int 
p{k)dk = 1 - {-iT\nt)-^+^ 

~ exp[-(-zTlnt)-'^+^], (11) 
an expansion valid for small t only. Given that (r|e~*^^) = ^(e"*''-' ), 

(rje--^) = exp[-(-.Tlnt)-/^+^], (12) 

at leading order as t 0. Substituting (fTT]) - (fT2|) into (ITOl) and making the change 
t ^ y = J2f=ii~j'^^^t) '^^^^ obtains: 



TV 



Y< = 1- exp 



-rina)-^+i5^j-^+i 



(13) 
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Figure 7. Numerically computed probability Y2 (averaged over 10'* disordered media) 
that two independent walkers choose the same site in the one dimensional problem with 
iV = 10^. The straight Hue is relation (PTjl. 



If < 2 (superabundant regime in Id [265), above sum diverges with A^: 



-'2 



1 — exp 

Thus, one may conclude that 

¥2 = 1 (1< /? < 2 



-f-T Ina) 



2-/?. 



1 (N^oo). 



(14) 



(15) 



hke in the deterministic case, for sufficiently large systems at any finite noise intensity T: 
The most attractive target at i* completely dominates the sum - and the probability 
of choosing i* among the N targets tends to unity even at finite T. 

On the other hand, in the disorder-free case (homogeneous medium or (3 ^ 00), 



E 



N 
1=1 



_2_ 



E 



N 
i=l 



1 

N 



0, 



(16) 



for any finite T > 0. Hence, all targets become practically equally attractive and each 
walker chooses one of them essentially at random (delocalized regime). 

For finite (3 > 2, numerical results indicate that Y2 remains finite as ^ 00, 
although it decays rapidly with (3, see figure [3 The values of Y2 plotted are practically 
unchanged for = 10^ and N = 10^. Hence, the finite values obtained for (3 > 2 can 
not be attributed to finite effects. 



6.2. Linear preference 

Another case of interest is the linear preferential choice f{x) = x, i.e., rj = ki/i. It is 
exactly soluble and closely analogous to the random energy model [36]. In the disorder- 
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free case {ki = 1), Y2 reads 
at large N. 

In the presence of disorder with (3 > 2, one has (rj) = — l)/(/5 — 2) < 00. 

Since only the small t behaviour contributes in the product in ffTOj) . then (exp(— trj)) ~ 
1 — t{ri) ~ exp(— t(rj)) and (r^^ exp(— tr,)) ~ (r,)^ exp(— t(rj)). Substituting these 
expressions in (|TOl) . one finds: 

/■oo / ^ \ --2 



the same result as the disorder-free case flTTl) . Therefore, disorder is irrelevant in 
heterogeneous media with (3 > 2. 

For 1 < /3 < 2, (rj) = 00. It is convenient to write 



/oo 
dk k-^[l-e^Y>{-tk/i)\. 



(19) 



Since dk k ^[1 — exp(— tfc/z)] is integrable and tends to as t — > 0, the lower bound 
1 can be replaced by in the integral of f|T9l) . One then obtains, 



(e-*^^-) = !-(/?- l){t/iY-^ j dx x-^{l - e"") 

^ exp [-(/3 - l)(tA)^-^(-r(l - /?))] , (20) 

for small enough t. By using (120|) and the identity (r|e~*''j) = ^(e"*''^), (fTOi) can be 
integrated and gives the simple result: 

F2 = 2-/? (l</?<2). (21) 

The results are confirmed numerically in figure [3 Like in the previous exponential 
preference, Y2 is no n- vanishing in the limit of large system sizes for strongly 
heterogeneous media {(3 < 2). This result is asymptotically identical to that obtained 
for the second moment of the weights of the microscopic states in the partition function 
of the random energy model at low temperatures [36]. It is also identical to the Y2 
parameter for a sum of positive Levy variables a problem equivalent to setting 
Tj = ki instead of = ki/i in (ITOl) . This shows that in strongly disordered environments, 
the distance travelled to the best target does not play a major role compared to the target 
size in the choice process. However, since 1^2 < 1, noise is not completely irrelevant: 
there is still a finite probability that the two independent walkers choose a different 
target. 
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7. Discussion and conclusions 

We have shown that a "knowledgeable" walker following simple deterministic rules 
in a medium composed of targets with heterogeneous weights describes trajectories 
characterized by multiple length scales. For a particular distribution of target weights, 
the sojourn lengths are maximally fluctuating and distributed in a Levy-like manner, 
as P{1) ~ with a ^ 2. As shown in [25], this model reproduces observations of 
spider monkeys foraging patterns, where a ~ 2 has been measured. These monkeys 
feed on trees whose size distribution is broad, approximately an inverse power-law with 
exponent f3 ~ 2.6 |25], not far from the special value /5c — 3 found in the model. 

The type of processes modeled here (but see also [35] for a similar approach) also 
represent a useful alternative to merely random-walk descriptions of human travels: 
humans use (mental) maps, take non-random decisions and concentrate their activity 
very heterogeneously in space. Qualitatively, our model could describe the movement 
patterns of hunter-gatherers [15] or the diffusion process of bank notes mediated by 
the travels of modern humans at the scale of a country [16]. The former system is 
empirically well described by a step length distribution with exponent a ~ 2 [15], and 
the latter by a ~ 1.6 p^. In |16j, long enough steps connect different cities and the 
power-law behaviour could be induced by the broad distribution of city sizes, that can be 
approximated by a Zipf's law [3Dj, i.e. relation ([T]) with j3 = 2. Even-though our results 
agree only qualitatively with these observations (no scaling behaviour in P{T) is actually 
found here for jS = 2), the present framework, appropriately modified, may provide a 
promising alternative for the description of this process. Future research should consider 
more sophisticated decision rules and systems of targets not uniformly distributed in 
space. 

On a theoretical point of view, the model proposed exhibits rich dynamical features 
in spite of its simplicity. We find, for example, that not every scale-free weight 
distribution induces walks with scale-free properties. Besides, the fact that the mean 
field approximation worked out in [26] is not very accurate, suggests that memory 
effects (arising the avoidance of previously visited sites) play a crucial role. The special 
exponent value a = 2 found in simulations is actually not well understood. 

We also present arguments showing that in the presence of strong disorder (small 
values of /3 here), the addition of stochasticity in the choice rules should not affect the 
main properties of the walks, provided that the probability of visiting a site increases 
sufficiently rapidly with the site attractiveness. This situation is similar to the lack of 
self-averaging in spin glasses at low temperature [37]. As shown by our calculation of 
Y2 in the one-dimensional case, noise should not affect much the distribution of sojourn 
length P{1), thus explaining a posteriori the results found in [25j in a variant of the 
model including forager mistakes. The analytical results presented here being limited 
to a single step process, it would be interesting to further investigate whether there is 
still a finite probability to find two (or n) independent walkers at the same site after 
many steps. Such probability does remain asymptotically finite in several random walk 
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models in disordered media, the remarkable Golosov localization phenomenon [S5] . 

A system of many non-interacting travelling agents with mental map are susceptible 
to make contact (when they occupy a same site) much more frequently than random 
walkers. The results above suggest that such contact networks could be, in some cases, 
controlled by the quenched randomness of the medium in spite of the stochasticity of 
the agent's decisions. Such robustness could have important social implications, for 
instance for the formation of social networks or the spreading of diseases [39] . 
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